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Abstract 

Dynamical Symmetry breaking and meson masses are studied in the Nambu-Jona-Lasinio (NJL) 
model at finite temperature and chemical potential using the dimensional regularization. Since 
the model is not renormalizable in four space-time dimensions, physical results and parameters 
depend on the regularization method. Following the imaginary time formalism, we introduce the 
temperature, T and the chemical potential, fi. The parameters of the model are fixed by calculating 
the pion mass and decay constant in the dimensional regularization at T = /u = 0. 
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I. INTRODUCTION 



QCD is a fundamental theory of quarks and gluons whose couphng constant is large 
at low energy scale. Therefore one cannot adopt the perturbative expansion in powers of 
the coupling constant at low energy. To study physics in a hadronic phase we can not 
avoid considering the non-perturbative effect of QCD. One of the possible ways to evaluate 
the phenomena in the hadronic phase is to use a low energy effective theory. Some phe- 
nomenological parameters are introduced to construct the low energy effective theory which 
is simpler than QCD to deal with. NJL model is one of the low energy effective theories of 
QCD [ij , for a review see Refs.QSflB. The model has the same chiral symmetry as QCD 
and the symmetry is broken down dynamically. Thus it is often used to study symmetry 
properties of QCD in the hadronic phase. 

The NJL model contains four-fermion interactions. Since a four-fermion interaction is a 
dimension six operator, the NJL model is not renormalizable in four space-time dimensions. 
To obtain finite expressions we must regularize the theory. In such a non-renormalizable 
model most of the physical quantities depend on the regularization method and a parameter 
introduced to regularize the theory. For example, a cut-off scale is introduced as a parameter 
in a cut-off regularization. The coupling constant and the cut-off scale are determined 
phenomenologically. 

The regularization in general may break some of the symmetries of the theory. For 



example a naive three dimensional cut-off breaks Lorentz and gauge invariance |6|]. In the 
present paper we employ the dimensional regularization. It preserves most of the symmetries 
of the theory, including the general covariance. We regard the space-time dimensions as one 
of the parameters in the effective theory. Thus the dimensions should be determined in some 
low energy phenomena which have some relation with chiral symmetry breaking. 

There are some works investigating NJL model in the dimensional regularization. The 



general properties of the renormalization and the renorma^ 



ization group is studied in arbi- 



trary dimensions 2 < D < 4 Iq]. In Refs.jlO, 11, 12, 13, Q,ll5|,Q the NJL model 

is considered as a prototype model of composite Higgs. The phase structure of dynamical 
symmetry breaking is analyzed at high temperature, density, electro-magnetic field and cur- 
vature. But we have no established fundamental theory of composite Higgs models at high 
energy scale. Thus the physical scale of the theory has not been fixed and the contribution 
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of the current fermion mass has not been analyzed. 



In Ref. 171] the dimensional regularization is modified to keep four dimensional properties 
of the non-renormalizable theory as much as possible. To achieve this goal the dimensional 
regularization is applied to only the radial part in loop integrals. It is one of the analytic 
regularization. The meson loop contribution to the chiral symmetry breaking is also analyzed 



in the NJL model with the modified dimensional regularization [18|. The procedure is easily 
extended to the finite temperature and chemical potential. However, the phenomenologically 
consistent dimensions are less than two in this approach. Because of the infrared divergence 
it seems to be difficult to obtain a finite result at finite temperature. 

In the present paper we regard the NJL model as a low energy effective theory of QCD 
and apply the dimensional regularization not only to the radial part but also to the angular 
parts of internal momenta in loop integrals. We include the effect of the current quark mass 
and fix the scale through the observed properties of pseudo-scalar mesons, the pion mass 
and the decay constant at zero density and temperature. The constituent quark mass and 
the meson properties are evaluated in thermal equilibrium. 

The paper is organized in the following way. First we apply the dimensional regular- 
ization to the NJL model with two fiavors of quarks. As is well-known, the small mass of 
quarks explicitly breaks the chiral symmetry. Thus the chiral symmetry is only approxi- 
mate. In such a model we evaluate the phase structure of the theory. We need to perform 
renormalization as well to define a positive coupling constant. In Sec. II we calculate the 
mass of scalar and pseudo-scalar mesons. We take the massless quarks limit and show that 
the dimensional regularization keeps the Nambu-Goldstone modes massless. A relationship 
between the space-time dimensions and the physical mass scale is also discussed. In Sec. IV 
we introduce the thermal effect in the imaginary time formalism and evaluate the effective 
potential. In Sec.V we calculate the mass of scalar and pseudo-scalar mesons at finite T 
and /i. The three massless poles of pseudo-scalar mesons survive at higher temperature and 
chemical potential. At the end concluding remarks are given. 

II. DYNAMICAL SYMMETRY BREAKING 

The chiral symmetry is broken when the composite operator of a quark and an anti-quark, 
ipip, develops a non-vanishing expectation value. It is caused by the QCD dynamics and is 
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called dynamical symmetry breaking. 

The NJL model is one of the simplest models of djTiamical chiral symmetry breaking. 
The model is originally introduced to evaluate the low energy properties of hadrons as a 
bound state of some primary fermion fields The Lagrangian of the two-flavor NJL model 
is defined by 

£ = ^(z^ - m)^ + gl {{i^i^f + (V^275r»'} , (1) 

where is an effective coupling constant, t"" represents the isospin Pauli matrices, and 
m = diag(m„,mrf) is the mass matrix of up and down quarks. In the Lagrangian ([T]) we 
omit the color and the flavor indeces. 

In the limit of massless quarks this Lagrangian is invariant under the global flavor trans- 
formation 

^ ^ e'e°">, (2) 

and the chiral transformation 

^ ^ e^^^-^^V. (3) 

The quark mass term violates these symmetries explicitly. When the composite operator 
iptp develops a non-vanishing expectation value, the quark acquires a mass and the chiral 
symmetry is broken dynamically. 

The vacuum expectation value of the composite operator ipip can be found by solving the 
gap equation. The gap equation of the NJL model in D dimensions is 

{a)=2^gl j ^^^S{kl (4) 

where S{k) is the fermion propagator and "tr" denotes trace with respect to flavor, color 
and spinor indices, 

S{k) = -. (5) 



Integration over k takes the gap equation to the form 

^2l(D/2-l) 

je{u,d} 

^ 9lA{D) J2 K + (^))[K + (6) 
je{u,d} 

where is the number of colors and A{D) is defined by 
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FIG. 1: Typical behavior of the effective potential {D = 2.8, niu = 3MeV, = 5MeV and 



c/° = -O.OlMeV^-^) 



The gap equation ([6]) has three solutions. To find the stable solution we evaluate the effective 



potential for the scalar channel, a, (see, for example Ref. 10|) 

(t2 _A{D 
4^~^ 



je{u,d} 

Extrema of the effective potential satisfy the gap equation The most stable solution is 
determined by the minimum of the effective potential. As is shown in the Fig{Tl we find the 
minimum at a negative value of a for a negative coupling constant, g^, with a finite current 

quark mass. We note that the expectation value, (ipip) ~ ni'^): has a negative value as 

usual. 

At the massless limit, m ^ 0, the gap equation ([HD is simplified and the solutions are: 

{a) = 0, (9) 

The non- vanishing solution ( ITOl) is stabler, see Fig.l. It shows that the composite operator 
i/jiIj develops a non-vanishing expectation value and the chiral symmetry is broken dynami- 
cally. It should be noted that 

A{D) < 0, for 2< D <A. (11) 

The left hand side of the Eq.l ITOl) is real and positive. To find a real solution of Eq.l ITOl) we 
must set a negative value for the coupling constant g^. It is one of characteristic features of 
the dimensional regularization. 



A. Renormalization 



The Lagrangian (1) is not renormalizable for D = 4. However, the renormahzation is 
sometimes useful to connect the results in the different regularization schemes. Here we 
renormalize parameters in our model. 

In the leading order of the expansion the radiative correction of the four-fermion 
coupling for the scalar channel is given by the summation of all bubble type diagrams, 



GsipMcr)) 



X*>0' 



+ 



+ 





0\2 



where the self-energy Uglp"^) is 

n.(/) = 




0\2 



(2^ 



tr[^(A;)^(A;-p)] 



(12) 



2{g:r{D - 1)A{D) J2 fd^L 



D/2-1 
3 ' 



(13) 



with Lj = {rrij + (cr))^ — x{l — x)p'^. 

We renormalize the coupling constant (7° by imposing the following renormalization con- 
dition 

Gs{p' = 0, (a) = Mo) = Zg{Mo)G:{p' = 0, (a) = Mq) 

4^.(Mo)^?° 



je{u,d} 



(14) 



j&{u,d} 

where we introduce the renormalization constant, Zg{Mo). The superscript r stands for the 
renormalized quantities and Mq is a renormalization scale. 



At the limit = the self-energy Uslp^) simplifies to 



Hsip' = 0)|(.)=M„ = 2{g:riD - l)AiD) ^ [(m, + Moff/'-'. (15) 

je{u,d} 



Thus the renormalization constant is given by 



(16) 



Therefore the renormalized coupling g^. is found to be 

1 1 1 



ie{«,d} 



+ {D-l)AiD). (17) 



We rewrite the gap equation ([6]) in terms of the renormalized coupling constant as 



Yl K- + (^))[K- + 



2]_D/2-l 



^"^"'^ MItt-T^ + ^-I)- (18) 



5^ [(m, +Mo)2]^/2-i \^9:AiD) 
je{u,d} 

In the massless limit m — 0, the non-trivial solution of the gap equation is given by 

((a)2)(^/2-i) 



( \— + D-l] (M2)(^/2-1)_ (ig) 

A real solution for (cr) exists for 1/(7^ < l/g^^, 

9n = \ ^ > for 2 < D < 4. (20) 

2(1-D)A(L)) ^ ' 

The critical coupling is real and positive. It should be noted that the broken phase is realized 
for gl > g'J>0 or g^ < 0. 

It may be useful to rewrite the Eqs. ( fTTj) and ( ITSl) in terms of the critical coupling, 

1 1 11 



je{u,d} 

and 

J] (m, + (a))[(m, + (a))T-^/^ 
j6{M} 2(a) / 1 1 



(21) 



J2 [{mj + Mof]'-""/^ MD) \g: g'J 

je{u,d} 



(22) 
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From Eq. (l2T|) it is clear that the bare couphng should be negative in the broken phase, 
1/(7^ < 1/(7^^, in the dimensional regularization. 

The renormalization group (3 function is defined by j^, 9], 



m) = Mo 



dMo 



Using Eq.(l22]) in Eq.(l23]) we obtain 

Mo [{m, + Mo)T-'^/' 



(5{gl) = 2{D - 2) 



j&{u,d} 



91' 



(23) 



(24) 



The ultraviolet stable fixed point appears at the critical coupling, = g'^ 



III. MESON MASSES 



To describe a bound state in the hadronic phase non-perturbative effects of QCD should 
be considered. Here we study meson masses in the NJL model at the leading order of the 
expansion. 



A. Pseudo-scalar mesons 



First we calculate the pion mass in the two-flavor NJL model. It corresponds to the 
Nambu-Goldstone modes at the massless quark limit. Since the current quark mass explicitly 
breaks the chiral symmetry, the pion mass should be proportional to the current quark mass 
in the broken phase. 

The pion mass is deflned by the pole of the propagator for a pseudo-scalar channel. The 
leading order of the expansion is described by the summation of inflnite number of the 
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bubble diagrams, 




+ 71= 



+ 7r= 





0\2 



cab 



2gl - nf (p2 



m^a — + 0{p 



(25) 



where we insert the renormahzation scale M^^^ ^ to impose the correct mass dimensions 
on Zt^. Jl^ij?^ is the self-energy for the pseudo-scalar channel which is given by 




tr[i7V"^(A;)i7V^(A;-p)]. 



(27r)^ 

Integrating over the D-dimensional momentum /c, we obtain 

+§(1 - x){^d + {a)f + (l - f ) + (cr))(m, + (a)) 
X [x(m„ + (a))^ + (1 - a;)(mrf + {a)f - x(l - x)p^] '^^^"^ 



(26) 



(27) 



and 



n^(p2) = 2(^°)^A(D) Trfx 5^ [{I - D)x{l - x)p' + {m, + {a)f] 
X [(m, + {a)f - x(l - x)p2]^/'-' . 



(28) 



We can analytically calculate the momentum integral at finite p^, but it is enough to 
evaluate the self-energy near = to observe the massless pole of the Nambu-Goldstone 
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mode. Straightforward calculations lead to 



n^''(/ ~ 0) = A{glfA{D) 



and 



2 y: [k + M)^/^-^ 

je{u,d} 



where is given by 



[(md+(a))2-(m„ + (a))2]3 
x{(|-l) [(m. + (a)r^-(m.+ (a))^+^] 



In the isospin symmetric case, rUu = ma, Ill''^{p'^ ~ 0) coincides with n^p^ ~ 0). 
Further we obtain the pion wave function renormalization constants, 



Z;,\ = -A{D)f{D)M, 



4-D 
' 



and 



je{u,d} 



Neglecting 0(p^) terms, the pion masses are found to be 



and 
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In the massless quark limit Eqs. (lM|) and (135|) read 

Z-'Mt'ml _ A{D){{a)r"-'. (36) 

Substituting the non-trivial solution of the gap equation (fTOj) into Eq. (!36|) . we observe that 
the pion mass disappears. Therefore the pion propagator has massless pole and thus the 
Nambu-Goldstone mode remains in the dimensional regularization. This is natural given 
the fact that dimensional re^larization does not damage Word-Takahashi identity, whereas 
cut-off regularization does [6]. If the constituent quark mass is not generated, the above 
expression can be written as 

W-V^4^(4-4)(^©--. (37) 

Since the renormalization constant is positive definite, the pion mass squared is negative 
in the broken phase, 1/(7^ < ^/ g'^ ■ Such a tachyon pole is compensated by the constituent 
quark mass. On the other hand inverse of the renormalization constant Z^^ is divergent 
and a massless pion pole appears in the symmetric phase, 1/(7^ > 1/ g"^ ■ 



B. The scalar meson 

Next we calculate the meson mass in the scalar channel. The propagator for the scalar 
channel in the leading order of the expansion is given in Eq. fll2l) 



The scalar meson mass (given that the corresponding pole is not located on the real axis) is 
obtained by observing the extremum for \Gs{p'^ = m"^, (cr))|. We can find the corresponding 
pole for D > 2.2. It can not be found for D < 2.2 and for soft mode in the next section. The 
extremum condition Eq. fl38l) for scalar meson mass is applicable for more general cases. 

First we discuss the massless pole in the scalar channel. It is not realistic but useful to 
understand the symmetry properties of NJL model. If the scalar meson is light enough, we 
can employ the same procedure as in the previous subsection. Thus we use the definition 
for rricr similar to m,r- 
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The self-energy Uslp^) defined in Eq. (fT3!l is expanded near ~ as 

n,(/ ~ 0) = 2{g:nD - 1)A{D) J] [(m, + (a))^/'^^ 



2]D/2-2 



O(p^). (40) 



Thus the renormalization constant for the scalar wave function is obtained by 

D 



^ (l - f ) A{D)Mt^ Y: [(-^- + (-))r^'-'- (41) 

^ ^ j&{n,d} 



It has a positive definite value for 2 < D < 4. In the present approximation the scalar mass 
is given by 

Z^'M^'ml = ^ - Ia{D){D -1) J2 [K- + (42) 

j€{u,d} 

In the massless quark limit it becomes 

Z^'M^'ml ^(2 - ^) = - ^) (2 - I^)(Mo^)^/-^ (43) 

where we used the gap equation (fTOl) . As is shown in Sec. II, the bare coupling is always 
negative in the broken phase for 2 < D < A. Thus the sigma mass squared is real and 
positive for l/g^. > ^/g^- The sigma mass decreases as the inverse of the coupling 1/g^ 
increases. For 1/g^ > l/g^ the renormalization constant Zcr vanishes. Hence a massless 
pole appears for the scalar channel. The sigma meson degenerates with the pion in the 
symmetric phase. 

If we do not consider the constituent quark mass, the sigma meson mass degenerates with 
the pion mass (|371) in the limit of massless quarks. The tachyon pole ( |371) shows that the 
vacuum is unstable and the constituent quark mass should be generated. 

Since the scalar meson is heavy in the real world, we do not apply the expression fH2|) 
to evaluate the real scalar meson mass. Furthermore, there is no pole for < 4(((t) + m)^. 
An imaginary part appears in the self-energy for p"^ > 4((cr) + m)^. The imaginary part 
corresponds to the decay width of the sigma meson. To find the sigma meson resonance we 
numerically calculate the pole of the propagator. 
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C. Dimensions and physical scale 



In this subsection we discuss the physical scale and coupling constant in terms of the 
space-time dimensions. Since the NJL model is phenomenological model of QCD, the scale 
of the model should be determined to reproduce the observable physical quantities. Here we 
calculate the pion decay constant and the pion mass as typical quantities in the hadronic 
phase. To consider the QCD interaction we set Nc = 3. In the present model the space- 
time dimensions D in the loop integrals is also regarded as one of the parameters of the 
phenomenological model. 

Inserting the measured values m° = 135.0MeV and = 139.6MeV into Eqs. fl34l) and 
f l55|) . we can obtain the matching condition for the coupling constant g,, and the space-time 
dimensions D with the physical scale. But the main part of the mass difference between 
the neutral and the charged pion arises from the electro-magnetic interaction which is not 
considered here [19i]. Thus we inspect only the isospin symmetric limit m° = ~ 136Mev 
and set rriu = below. 

We start with the relationship in the current algebra (See, for example, appendix in Ref. 
{2^.) The pion decay constant is written as 

h = -Z-'''{a). (44) 



From the Gell-Mann-Oakes-Renner relation 



24| we obtain 



/>^ = ^Mo^-^(m„ + m,). (45) 

The constituent quark mass (cr) is determined by the gap equation (jS]). Substituting the 
measured values = 136MeV, /,r = 93MeV and the solution of the gap equation we obtain 
the matching condition for the coupling constant g^^, the space-time dimensions D and the 
renormalization scale Mq with the physical scale. After some numerical calculations we get 
the physical scale for (7,^ and Mq as a function of the space-time dimensions D. 

In Fig. 2 the curves satisfy the relations (jBl) and ( l45i) o\i D — g^ and D — Mq plane for 
rriu = TTLd = 3, 4 and 5MeV. As is shown in Fig. 2, the bare coupling constant and the 
renormalization scale has no strong dependence on the current quark mass. The renormal- 
ization scale is found to be below the pion mass. In Fig. 3 we illustrate the behavior of the 
constituent quark mass and of the expect value of the composite operator of a quark and 
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mu=m(j=3MeV 
mu=mci=4MeV 

mu=m(j=5MeV 



> 

CD 



mu=m(j=4l\/leV 




2 2.5 3 3.5 4 
D 

FIG. 2: The physical scale for the bare coupling constant and the renormalization scale Mq are 
shown as a function of the space-time dimensions D for Nc = 3, rriu = rrid = 3, 4 and 5MeV. 
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2 2.5 3 3.5 4 
D 

FIG. 3: The constituent quark mass (cr) and the expectation value {tpil^) are shown as a function 
of the space-time dimensions D. 



an anti-quark, {ipip) 



(a), for the parameters fixed by and fHSl) . To obtain the 



constituent quark mass near 300MeV we should consider the lower space-time dimensions 
or larger current quark mass. The expectation value {ipip) has strong dependence on the 
space-time dimensions. But the normalized value, {ijj^p)MQ~^ , varies less than 0.1% between 
two and four dimensions and is found to be 



(mm)M( 



4-D 



-(298.7MeV)3 (m„ = ma = 3MeV), 
-(271.4MeV)3 (m„ = = 4MeV), 
-(252.0MeV)3 (m„ = ma = 5MeV). 



(46) 



Inserting the non-trivial solution of the gap equation into Eq.f l35l) and taking = 
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mu=mci=5MeV 
mu=m(j=4MeV 
mu=mtj=3MeV 



mu=md=3MeV 
mL,=mcj=4MeV 
mu=mcj=5MeV 
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D 



3.5 




756 758 760 762 764 
p (MeV) 



(a) Pion mass from Eq.jSS]) (dashed (b) G^^{p^, {a)) for D = 2.4 and 
line) and Eq.(|48]) (full line). rriu = rud = 5MeV. 

FIG. 4: Solutions for Eq.l^ and Eq.diOD. 
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(a) Sigma meson mass from Eq. (09]) (b) Sigma meson width from Eq. §1 
(full hue) and Eq. ([50]) (dotted hue). 

FIG. 5: Sigma meson mass and width 



m we get 



m 



+ 



m 



(47) 



The Gell-Mann-Oakes-Renner relation fH5l) is reproduced from Eqs. fH^j) and fl47l) . Thus the 
Eq. (135!) should be consistent with the model independent relationships (jH]) and (H5l) . 

In Figs. 4 and 5 [4^ we calculate the pion and the sigma meson mass by numerically 
solving the equations 

1 1 



G^\p' = mi, {a)) 



2^?° 4(^?o)2 



0, 



(48) 
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and 

G:\p' = ml (a)) = ^ - [Usip' = ml] = 0, (49) 

using the gap equation. As is clearly seen in Fig. 4 (a), the results for are found at the 
physical pion mass scale ~ 136MeV for D > 2.2. The approximation in Eq. ( |35ll seems to 
be valid for higher dimensions. A typical behavior of the real and the imaginary part of 
Gj^Ip"^, (cr)) is shown in Fig. 4 (b). Since the sigma meson is almost stable in the leading 
order of the expansion. We use the approximate equation (H9|) . The approximation 
is not valid in some cases, for example, D < 2.3, hight temperature and/or large chemical 
potential. In such cases we use 

d_ 

dp' 

The sigma meson develops the heavier mass for higher dimensions in Fig. 5 (a). If imaginary 
part of Gs is small, we can identify the sigma meson width, Fo-, as 

r. = -|^^im[n3(/ = m^)]. (51) 

with TTio- obtained from Eq. We plot the behavior of the sigma meson width in Fig. 5 

(b). For smaller dimensions both the sigma meson and width approach to the result obtained 
by solving the Roy equations in 2l(]. As is pointed out in Ref 22(], higher order corrections, 
TTTT scattering, are essential for the sigma width. To discuss the sigma meson width one 
should calculate the next to leading order corrections in expansion. 



\G,ip' = ml{a))\ = 0. (50) 



IV. DYNAMICAL SYMMETRY BREAKING AT FINITE T AND // 

To define the thermal equilibrium the time direction has to be fixed. This breaks general 
covariance however, it is important not to introduce an additional violation of the covariance 
through the regularization. Thus we apply the dimensional regularization to the NJL model 
in the thermal equilibrium and extend the above analysis to the finite temperature T and 
the chemical potential /i. 

First we discuss the phase structure of the theory. It is expected that the broken chiral 
symmetry is restored at high temperature and/or large chemical potential. Following the 
standard procedure of the imaginary-time formalism, we introduce the temperature and the 
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chemical potential to the NJL model in D dimensions 10|. Hence the gap equation (jlj) is 
modified as 

(52) 



where Spf^i^k) is the fermion propagator at finite T and /i, 

1 



(53) 



■ 7 - (a;„ - ■j/i)74 + m + {a) - ie 
Due to the anti-periodicity of the fermion field the A;° integral is replaced by the summation 

over Un = {2n + l)7r//?. We note that the fiavor symmetry is broken by the quark mass and 

the chemical potential in the thermal equilibrium (see, for example, Ref. 

Integrating over the momentum in Eq. fl52l) . we obtain 

(a) = -V2g:A{D - 1) Yl + (^)) 



23|.) 



,21 (D-3)/2 



n 

We also perform the summation in Eq. (l52|) and get another expression, 

TJij + {a) 



(54) 



je{u,d} 



(27r)(^-i)/2^ fD 



1 /"oo . 



X 



1 



l + e 



/3(^fc2+(m, + (<7»2+^) ^ _^ g/3(^fc2 + (m,+(<7»2-|.) 



(55) 



10| at the massless quark limit. 



The result exactly reproduces the one obtained in Ref. 

At the zero temperature limit the gap equation (155!) simplifies to 

rrij + {a) 



je{u,d} 



(2vr)(^-i)/2"^ /D 



I- 1\ ^ 



iG{«,rf} 



xe{fi - ^Jk^ + (m, + (56) 

For < m„ + (cr) and fi < rrid + (cr) the theta function vanishes and we obtain the same 
solution with zero chemical potential. 

je{u,d} 
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(a) Behavior of the effective potential (b) Behavior of the constituent quark 
at ;U = 300(MeV). mass (a). 

FIG. 6: Typical behavior of the effective potential and the constituent quark mass at finite T and 
fi in D = 2.4 and rriu = rrid = 5MeV. 

Hence, the chemical potential does not contribute to the constituent quark mass (a) at 
T = 0, yU < rriu + (o") and /i < + (cr). 

To find the stable solution we evaluate the effective potential V{a) at finite T and /i. It 
is obtained by integrating over (a) in the gap equation fl55l) with further replacement (a) 
by a 



V{a) 



A{D) 

1 y/2N^ 1 
~^(27r)(^-i)/2 fD-l 



dk 

j(^{u,d} 



(58) 



X I In , h In 



The stable solution of the gap equation is obtained by observing the minimum of the effective 
potential. We numerically calculate it and show typical behavior of the effective potential 
V{{cr)) and of the dynamically generated quark mass (a) in FiglHl A first order transition 
takes place for n = 300 (MeV). In the massless quark limit a phase transition occurs as 
is shown in Ref . jiol] . Because of the explicit chiral symmetry breaking term, i.e. finite 
current quark mass, the phase transition becomes a cross over in Figj6] 25, 26|. In Ref. [lo| 
detailed behavior of the dynamically generated mass and the phase structure are found for 
TTLu = fUd = 0. Analytic expressions for some characteristic points are also presented there. 
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V. MESON MASSES AT FINITE T AND fi 



Next we consider meson masses in a thermal equilibrium. The question of whether the 
Nambu-Goldstone modes of the chiral symmetry breaking are still massless or not at finite 



T and yU is not trivial 



27 



28 



30, l3l| • It is expected that the dispersion law is modified 



in the thermal equilibrium. 

In the imaginary time formalism the self-energy for the pseudo-scalar channel, Il'^{p'^) is 
modified as 



Performing the trace operation in Eq. fl59l) . we obtain 



(27r)^-i 

{un - ifi){uJn - P4 - ifi) + k ■ {k — p) + {niu + {(T)){md + ((x)) 



+ {u ^ d), 



and 



(60) 



(61) 



je{u,d} 



{un - i^J-){uJn - Pa - ifi) + k ■ {k — p) + {rrij + {a) f 



[{un - ifJ^y + Elk ] [(t^n -Pa- if^y + E. 



2k 



where 



E- 



Ik 



e + {nij + {a)f, E2k^ = {k- pf + {m.j + {a)f. 



(62) 



We observe the massless pole of the pion propagator at the massless quark limit, again. For 
'f^u = = Eqs.(|60ll and (|6Ti) reduce to 



d^-^k 



(2n 



X 



1 



-inY + El 



m=0 



-Pa- ifJ'Y + E2 



m=0 



pI + p^ 



[{uJn - iyU)2 + Ei\m=0 ][{uJn - Pa - ^/i)^ + -^2 



m=0 



(63) 



To find the massless pole we put p"^ = p\ + p^ = 0. Then the final line in Eq. (!63!) disappears. 
Since the pseudo-scalar field with the momentum p obeys the Bose-Einstein statistics, the 



19 



fourth element is written as = 217t/(3 {I is an integer). We can simply change variables 
Un — Pa and k — p in the second line of Eq. fl63l) to u;„ and k respectively. 
The pion self-energy (jHSD reduces to 



n^(p^ = 0) = 8iV.(^°)^2^/^lj] 



(a) /3 



(27r 



(64) 



The pion mass at finite T and /i is found to be 



z-'mi{p^ = mt' = 7^ - 



1 



n-(p2 = 0) 



'^91 4(^^0)2 

1 1 1 

2^+(^^4- 



(2vr; 



-trS^^(A;). 



(65) 



Substituting the gap equation (152|) into Eq. (l65!) . we find that the pion is also massless at finite 
T and yU. The pion is still massless if the quark state is represented by a non-trivial solution 
of the gap equation. It is noteworthy that we do not use the Feynman parametrization, it 
is not applicable to calculate loops at finite T and /x. 

Next we evaluate the contribution of the finite current quark mass. Here we consider the 
isospin symmetric limit for simplicity and set = = m again. In this case a pole of 
the pion propagator is given by 







1 n°(p4,P 
1 



1 



4(^?°)^ 



(cr) + m 



13^ 



d^-^k 



X ■ 



Pt + P 



(66) 



[{un - + ^i'][(^„ - P4 - ^^^)' + E2^] ' 
where we have used the gap equation ( |52l) in passing from the first line to the second and 
third one. Performing the summation in Eq.([HS]) 231], we get 







1 

2^ 



m 



iV,2^/2(p2+p2 



k 1 



n_(E2 



4 J (27r)^-iEi^2 

n_(Ei) -n_(E2) 



ipi - El - E2 ipi + El - E2 

n+{Ei) - n+{E2) I - n.{Ei) - n+^E^) 



ip4^ - Ei + E2 



ipi + Ei + E2 



(67) 
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where 



n±{E) = 



(68) 



To find the dispersion law in the Minkowski space we perform Wick rotation, ip^ —>■ po 
Thus the Eg. (1671) reads 







1 

2^ 



1 



(cr) + m 



N,2^I\pI-P^) f d^-^k 1 



4 



(27r; 

1 - n+{Ei) - n_{E2) _ n.jEi) - n.{E2) 
Po- El- E2 po + El- E2 

n+{Ei) - n+{E2) 1 - - ^+(^2) 



D-l 



E1E2 



(69) 



Po- Ei + E2 po + Ei + E2 

If we change variable to p — fc in Eq. (l69ll . Ei and E2 are exchanged. Thus Eq.(l69 
simplifies to 







X 

+ 



1 



{a) + m 



+ 



N,2^I\pI-P^ 



{271 



,D-1 



^^^^ I - n+{Ei) - n^{Ei) 



El E2J pl-{Ei + E2y 
-n+iEi) - n4Ei) 



1 



1 



El E2J pI-{Ei-E2Y 
We divide the momentum integral in the radial and the angle parts, 



/•OO /•TT r'K 

/ k^-'^dk sin^-^ OidOi ■ ■ ■ sine 
Jo Jo Jo 



D~3d9D-3 



1-K 



de 



D-2- 



(70) 



(71) 



We can choose the space coordinate in the direction p = (p, 0, 0, ■ ■ ■ ), then k p = kpcosOi. 
The integrand in Eq. fl69|) depends only on k and 9i. The angle integrals with respect to 
6*2 ■■ ■ Od-2 give 



sin^-^ e2de2 ■ 

Thus the Eq. (|70l) reduces to 

1 







-\ 



sin O^^^dOij^ 



2n 



de 



-1 



D-2 



T{D/2-l)' 



(72) 



{a) + m 
2N,{pI-p') 



(27r)^/2r(Z}/2- 1) 



k'^-^dk I sin^-^^id^ 



D-3, 







X 



^^^^ I - n+{Ei) - n.iEi) 



El E2J pl-{Ei + E2Y-ie 
-n^Ei) - n^jEi) 
El E2J pl-{Ei-E2)^-ie^ 



^ 1 1 



(73) 
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FIG. 7: Dispersion law for pion at /i = for D = 2.4 and = md = 5MeV. 
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(a) Full line: Pion mass, dotted line: (b) IG5I at ^= lOOMeV. 

-2({a)+m). 

FIG. 8: Typical behavior of the pion pole and \G^\ at finite T and /i for D = 2.4 and niu = ma = 
5MeV. 

In the chiral limit m 0, the solution of this equation is Po — p"^ = 0, as is shown before. 
A finite current quark mass changes the solution. 

In the limit fi —>■ and T — > 0, the right hand side of Eq. (l73l) is a function of p"^. To see if 
the finite T and fi modify this dispersion law we numerically solve Eq. (!73l) . In the numerical 
analysis we put D = 2.4 and use the parameters determined in the subsection III.C. In Fig. [7] 
we plot the solutions on po,p plane at finite T. The solutions are on hyperbolic curves as in 
the usual T = case. The value of po on the curve at p = corresponds to the pion mass. 
The pion mass is obtained by solving Eq. (I73|) after rotating the contour by —tt/4. A soft 
mode for the pion channel is found by observing an extremum of the absolute value of the 
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propagator | Gs | . The pion mass decreases as T increases from to the critical temperature 
Tc. A soft mode appears above Tc. We also draw the behavior of IG5I near Tc in Fig. [8] (b). 
The sharp peak structure is observed near —2{{a) + m) for T > T^. In Fig. [S] (a) we see a 
gap at T ~ 60(MeV) for fi = 300(MeV). The gap comes from the one for (a) in Fig. [H] (b). 
For the scalar channel the self-energy, Ils{p4,p) is given by 

n.(p4,p) = / £^MSik),,Sik-p),,]. (74) 

We perform the trace operation and get 

n.fe, p) = 4iv.M)'2-/'i y: J £^ y: (75) 



X 



je{u,d} 

- ift){ujn - Pi - in) + k ■ {k — p) - {rrij + (cr))^ 



It is related to the self-energy for the pseudo-scalar channel, 



n.(P4,p) = n^(p4,p) 

- Y: 8((a)+m,)^iV.(^°)^2^/^^$: / ^"""^ 



(2vr) 



D-l 



je{u,d} 

[{un - ifiy + Eik"^] [{un -Pa- if^y + ^2fc^] ^"^^^ 

The difference between 115(^4, p) and 115(^4, p) disappears at the chiral limit ((a) +mj) 0. 
If the chiral symmetry is restored, the scalar and the pseudo-scalar channel degenerate even 
for finite T and fi. 

The two-point function for the scalar meson can be calculated similarly to Eq. (63) for 
the pion pole. After the summation over the Matsubara frequencies and the Wick rotation 
we obtain 



(a) + 

N,2^/''{pl-p^-A{{a) + m)^) f d^'-^k 1 



4 J (271)^-1^1^2 

I - n+{Ei) - n_{E2) n^{Ei) ~ n.^E^) 
Po- El- E2 po + El- E2 

n+{Ei) - n+{E2) 1 - n_(Ei) - ^+(^2) 
P0-E1 + E2 P0 + E1 + E2 



(77) 
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(a) Full line: sigma meson and soft 
mode masses, dotted line: 
-2{{a)+m). 

FIG. 9: Typical behavior of the sigma and pion pole at finite T and for D = 2.4 and rriu = ma = 
5MeV. 

Adopting the mathematical trick used in Eg. ( 1701) and integrating over the angle variables, 
6*2 ■■ ■ Od-2, we get 



(a) + 
- 4((o-) + m)^) 



(27r)«/2r(D/2 - 1) 







sm 



D-3 



+ 



_1 n l-n+(Ei)-n-(Ei) 
El E2) pl-{Ei + E2Y-ie 
-n+{Ei) - n_{Ei) 



It 



(78) 



Ex E2J pI- {El - r,2) 

The sigma meson and soft mode masses are obtained by observing the extremum of 
\Gs\ after rotating the contour by — 7r/4. We draw the behavior of the sigma meson mass 
and the soft mode in Fig. [9] (a). As temperature increases from to Tc, the sigma meson 
mass decreases for /i < 200(MeV). The sigma meson mass turns into the soft mode at 
temperature above T^. A gap at yU = 300(MeV) and T ~ 60(MeV) in Fig. [9] (a) is induced 
by the gap found in the < a > behavior shown in Fig. [6] (b). The other small gap appears 
at T ^ lOO(MeV). It comes from the extremum of < cr > as is shown in Fig 6. (b). In Fig. [9] 
(a) we also draw — 2((cr) + m) where a naive threshold (at double constituent quark mass) 
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for ipip state appears. The sigma channel above Tc we call the soft mode. In Fig. 9 (b) we 
show the extrema for \Gs\ which corresponds to the soft mode. 

The sigma and pion mass has been analysed in 3^ at finite T and /i in the cut-off 
regularization. For small /i a similar behavior is observed across the crossover transition. 
For larger /i the first order phase transition takes place in both regularizations. In this case 
the constituent quark mass (cr) shows some different behavior. For instance, the small peak 
in the constituent quark mass appears at yU = 300 (Me V) and T ^ lOO(MeV) in Fig. [9] (a) in 
the dimensional regularization. 



VI. CONCLUSION 



We studied characteristic features of the NJL model in the dimensional regularization. 
The dimensional regularization is applied to momentum integrals for internal fermion lines 
as usual. Since the model is not renormalizable, we can not take the four dimensional limit. 
We evaluated some physical properties of the model in the space-time dimensions less than 
four. We take notice of that only the radiative corrections should be evaluated in the space- 
time dimensions less than four to keep the four- dimensional properties in the real world. 
The classical parts are evaluated in four dimensions. 

In the NJL model we see that the approximate chiral symmetry is dynamically broken for 
a negative bare coupling. The constituent quark mass is also negative. After the renormal- 
ization, we can define the renormalized coupling constant. The dimensional regularization 
keeps most of symmetries. We have inspected the massless pole of the Nambu-Goldstone 
mode of the chiral symmetry breaking at the massless quark limit. Calculating the meson 
masses in the vacuum, we show that three pseudo-scalar modes have massless pole in the 
broken phase. It is the direct consequence of the Ward-Takahashi identity. In the dimen- 
sional regularization the Ward-Takahashi identity remains usual. We also see that the mass 
for the scalar mode is real and positive. 

The coupling constant, the space-time dimensions and the renormalization scale can be 
fixed by the pion mass and its decay constant. The dimensions less than three, but more than 
two, are preferable to obtain the constituent quark mass about 300 ~ 400MeV. The current 
quark mass has a non- negligible dependence on it. If we take the larger current quark mass, 
the constituent quark mass about 300MeV is realized for higher dimensions. However, a large 
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chemical potential induces the first order phase transition for lower dimensions 2 < D < 3 in 
;he massless quark limit, as is the case in other approaches, such as the cut-off regularization 



Schwinger- Dyson equation 
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36| and the lattice QCD [37 1. 



After we fix the parameters phenomenologically, we numerically evaluate the constituent 
quark mass at finite T and fi. The first order transition is observed for D = 2.4 at /i = 
300MeV. We also evaluate the thermal influence on the meson masses in the leading order of 
expansion. We show the dispersion law for pion and the behavior of meson masses. For 
lower temperature only a small influence is observed on the meson mass. It should be noted 
that boson loops have O(T^) contribution at finite temperature. A meson loop appears in 
the next to leading order of expansion. Though it is neglected in the present paper, 
the contribution from a meson loop may not be negligible at higher temperature. 

Some problems remain. One should include the QED corrections to evaluate the 
neutral and the charged pion mass difference jsS]. The color superconductivity is also in- 
teresting to consider in the extended NJL model by using the dimensional regularization 39|. 
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